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the presence of the nontrivial Haar measure in the SU(3) case, which has a stabilizing effect 
on the complexified dynamics. The freedom to modify and stabilize the complex Langevin 
process is discussed in some detail. 
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1 Introduction 

Because the fermion determinant in QCD at nonzero baryon chemical potential is com- 
plex, standard lattice QCD algorithms based on importance sampling cannot be used. As 
a result, a nonperturbative study of the QCD phase structure in the temperature - chem- 
ical potential plane is still missing [1]. Complex Langevin dynamics [2, 3] may provide a 
solution, since it is not based on importance sampling, but instead on a stochastic explo- 
ration of an enlarged (complexified) field space. However, the method is not guaranteed 
to work: the numerical solution of a complex Langevin process may converge to a wrong 
answer. This problem was observed immediately [4, 5] after complex Langevin dynamics 
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was proposed in the early 1980s and is still present in current complex Langevin simu- 
lations [6-8]. However, recently it has also been shown convincingly that in some cases 
complex Langevin simulations converge to the seemingly correct answer, even when the 
sign problem is severe, i.e. in the thermodynamic limit of one-, three- and four-dimensional 
field theories at nonzero chemical potential [9-12]. Given the importance of understand- 
ing strongly-coupled theories such as QCD at nonzero chemical potential numerically, this 
issue clearly needs to be addressed. 

The outstanding questions are therefore (1) to quantify whether the results from com- 
plex Langevin simulations have converged correctly; (2) to understand why the (in)correct 
convergence occurs; (3) to find a cure in the case of incorrect convergence. Recently we 
have clarified question (1) in some detail by putting the mathematical foundation of com- 
plex Langevin dynamics on firmer footing and deriving a set of criteria for correctness that 
need to be satisfied [13, 14]. These consistency conditions can be calculated during the 
complex Langevin simulation. 

In this paper we focus on question (2), yielding insight that can be used to address the 
third question. This is done in the context of two three-dimensional spin models recently 
studied: the abelian XY model and the nonabelian SU(3) spin model, both at nonzero 
chemical potential. The XY model was studied in ref. [8]. Here it was shown that correct 
convergence is obtained deep in the ordered phase, but incorrect convergence was found in 
the disordered phase and the transition region. This conclusion was based on a study of 
the expected analyticity of observables around fi 2 = 0, of properties of distributions in the 
complexified field space, and from a comparison with an alternative world-line formulation, 
which is sign-problem free [15]. The SU(3) spin model was studied in ref. [12], extending 
the classic papers [16, 17]. Using similar criteria as in the XY model (analyticity and 
Taylor series expansion around fi 2 = 0, localized distributions in the complexified field 
space) as well as a test of the criteria for correctness mentioned above, we concluded that 
complex Langevin dynamics yields the correct results in the entire phase diagram. This 
conclusion was subsequently supported by a study using a dual formulation, which is again 
sign-problem free [18, 19]. 

In the present paper we address this observed difference in performance of the complex 
Langevin method between the XY and the SU(3) spin model. Surprisingly, we find that 
the nonabelian nature of the SU(3) spin model is crucial. In particular we demonstrate 
that in the disordered phase of the XY model the real manifold is unstable against complex 
fluctuations while in the SU(3) spin model it is stable. Stability of the real manifold is 
important for understanding the expected analyticity of observables around fi 2 = 0. 

The paper is organized as follows. In sec. 2 we remind the reader of the SU(3) and 
XY models and reduce them to effective one-link models with complex couplings. In sec. 3 
the effective one-link models are discussed in detail and the difference between the abelian 
and nonabelian models is stressed. In sec. 4 we use the insight to illustrate how coordinate 
transformations may stabilize the dynamics. Conclusions and an outlook are given in 
ref. 5. Appendix A contains a detailed discussion of possible modifications of the complex 
Langevin approach. 
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2 SU(3) and U(l) spin models 



2.1 SU(3) spin model 

We consider the three-dimensional SU(3) spin model at nonzero chemical potential, with 
the action 

S = S B + S F , (2.1) 

where 

3 

s b = -p E E ( Tr u * Tv U U> + Tr u ^ u *+>) ' ( 2 - 2 ) 

X v = \ 

S F = -hY^ (e M Tr U x + e^Tr C/j) . (2.3) 

X 

The matrices U x are elements of SU(3). The action is complex and satisfies the usual sym- 
metry, S*(n) = S(—/JL*). The 'fermion' contribution follows from the full QCD determinant 
by preserving only propagation in the temporal direction at leading order in the hopping 
expansion (i.e. heavy quarks). The model is part of a whole family of effective Polyakov 
loop models with heavy quarks [20, 21]. A detailed complex Langevin study can be found 
in ref. [12]. 

We want to construct a one-link model that captures the essential dynamics. We 
therefore focus on a site x and consider the interaction of U x with its nearest neighbours. 
Treating all six neighbours equally and denoting their contribution as u € C, we write the 
effective one-link model as 1 

S u = -p 1 TrU- foTrU- 1 , (2.4) 

where the effect of the neighbours is captured, to a certain extent, by a simplified parametriza- 
tion 

ft = Q/3u + he", f3 2 = 6f3u* + he~", (2.5) 

with ft(n) = foi-fj*). 

Typical values of ft^ are determined by (3, h and \i, and by the contribution from the 
nearest neighbours, represented by u, in the original three-dimensional theory. We first 
note that the critical (5 value in the SU(3) spin model at vanishing fj, is /3 C ~ 0.133 [12, 19]. 2 
We therefore consider j3 values around this value, as in refs. [12, 19]. Writing u = \u\e l ' y , the 
correlation between its phase and amplitude can be determined using simulations in the 
full, three-dimensional theory. The result is shown in Fig. 1, from a study on a 12 3 lattice 
for 15 combinations of (3 and fj, at fixed h = 0.02. Note that during a Langevin simulation, 
u takes values in the familiar triangular shape with corners at 3e 2ni i/' 3 (q = 0,1,2). This 
explains the three spikes at ^y/n = 0, ±2/3. At larger f3 and/or fx, u lies predominantly 
in the trivial direction {q = 0) and 7 is closer to zero. Note also that at nonzero /i, 

x As always, U^ 1 is written instead of U* to allow for the correct extension to SL(3,C). 
2 Its relation with the four-dimensional coupling can be ultimately be understood from a combined 
strong-coupling/hopping parameter expansion [20]. 
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Figure 1. Correlation between the phase 7 and the modulus 6/3|u|, where u = \u\e 17 = TrU 
in the three-dimensional SU(3) spin model, for 15 combinations of j3 = 0.125,0.13,0.135, and 
H = 0.5, 1, 2, 3, 4, at h = 0.02 on a 12 3 lattice. 

the dynamics takes place slightly outside SU(3) [12] and \u\ is not strictly bounded by 3. 
Typical values of h and \x are determined by the relation with four-dimensional theory [22] : 
for Wilson quarks, h = {2k) Nt with n ~ 0.12 and N T = 4,6,8, etc. Relevant values of h 
are therefore small. Upon identification with the four-dimensional theory, \i as used here 
in fact corresponds to fx/T (with 1/T = a T N T ). 
We finally express the complex couplings as 



Px = /W 7 + he", h = {3 eS e-^ + he^, 

where /3 e g = 6/3\u\ < 2.5. 

2.2 XY model 

The action for the three-dimensional XY model at nonzero chemical potential reads 



(2.6) 



S = -f E E { etlS "' 3U * U Ui> + e~^UlU x+i ) 

X u=l 

= -p ^ E cos ^ x ~ ^ x+i/ ~ ^^3) , 



(2.7) 



where U x = e %l ? x are in this case U(l) phase variables. A complex Langevin study of this 
model can be found in ref. [8]. 

We first note that deep in the ordered phase, where <f> x — 4> x +c <C 1, the action reduces 

to 

^^^(^-^f, (2.8) 
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with /3x,2 = /3, /?3 = /3cosh/x. Here we dropped an overall constant and employed the 
periodicity of the lattice. Since this action is real, the sign problem is absent in this limit. 
Indeed, in the corresponding part of the phase diagram (large (5 and/or //), no problems 
were encountered in ref. [8]. 

As in the SU(3) case, we consider an effective one-link model of the form 

S u = -/3 1 U-(3 2 U-\ U = e i t (2.9) 

with complex couplings /Si 2- If we follow the approach from the SU(3) model and replace 
the six nearest neighbours by a common contribution u £ C, we find that 

Pi = (3u* (2 + cosh n) , fa = /3-u(2 + cosh / u), (2.10) 

so that /3j" = /?2 and the action Sjj is in fact real. Hence we are forced to represent the 
nearest neighbours by independent complex phases and we consider the general effective 
U(l) one-link model, 

S v = -faU-faU' 1 = -/3 lC os0-/3 2 sin0, (2.11) 
where we take P\ 2^P\2 ^ ^ without further restrictions. 

3 Complex Langevin for effective models 

3.1 U(l) one-link model 

We consider complex Langevin dynamics for the one-link models constructed above. In 
the U(l) case, the partition function reads 

Zu(i)= r d4e~ s ^\ (3.1) 

J — 7T 

and the corresponding complex Langevin equation is 

j> = K[<f>(t)]+v(t), K ^ = ~^ ^ 

with standard relations for the noise, (rj(t)) = 0, (Tj(t)rj(t')) = 25 (t — t'). Since the action 
and the drift term K are complex, <fi will take values in the complex plane and U = e l< ^ is 
no longer a phase variable. 

3.2 SU(iV) one-link models 

In the case of the SU(3) [or more generally an SU(iV)] spin model, we can express all 
dynamics in terms of the eigenvalues of U, subject to the constraints, i.e., 

= diag(e i * 1 ,e i * 2 ,...,e i M , <f>i + <f) 2 + . . . + 4>N = 0. (3.3) 



The partition function then includes integrating over the reduced Haar measure, which 
represents explicitly the integration over the full group manifold, 

a<b \ 1 J 
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(a, b = 1, ... , N) and reads 

^su(tv) = [ W #1 • • • d(j> N S (0i + 2 + . . . + 0jv) fl-({0 o }) e -^«*»» 

«/ — 7T 

/7T 
#1 . . . dcf> N S (0i + 2 + . . . + 0jv) e-^'M) , (3.5) 
-7T 

with 

S c s = Su + S H = Su - In F. (3.6) 

The non-abelian Haar measure, which is explicitly seen here as the reduced Haar measure, 
is the crucial difference between SU(iV) and U(l) spin models. This is especially clear in 
the case of SU(2) with a simple action 

Su = -^TrU = -/3cos0, (3 e C. (3.7) 

As in the U(l) case, the partition function involves a single integral over = 0i = —02 
only and reads 

Z S v(2)= rd0sin 2 (0)e-^W. (3 . 8) 

J — 7T 

The presence of the reduced Haar measure is the important difference with Eq. (3.1). 

To write down the Langevin equations in the SU(iV) case, we may follow several 
routes. Firstly, we may eliminate 0tv using the constraint in Eq. (3.3), and write Langevin 
equations for the remaining N — 1 degrees of freedom, 

a = K a [<j>(t)} + Va (t), K a {4>) = -^f-. (3.9) 

Each update requires N — 1 stochastic kicks. This is the approach followed in refs. [12, 16]. 

Secondly, the constraint can also be implemented by introducing new variables z a 
(a = 1, . . . , N), according to 

0a = z a - — (zi + z 2 + ■ ■ ■ + z N ) , (3.10) 

such that the constraint ^ a a = is automatically satisfied. This is a singular transfor- 
mation; the validity of the procedure therefore requires some justification, which is given 
in Appendix A. 4. 3. It should be noted that the stochastic process, even when restricted to 
the first N — 1 variables a , is different from the one above. We then write the following 
dynamics for z a , 

Z a = K a [ Z {t)\ + Tfe, K a (z) = (3.11) 

Oz a 

where the a are always considered as functions of the z a , according to Eq. (3.10). In this 
formulation N stochastic variables are used. However, it is easy to see that the real part 
of the centre-of-mass coordinate is freely diffusing (the imaginary part is not evolving at 
all), since 
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Since the real parts are angular degrees of freedom, this does not cause harm in terms of 
runaways. 

Finally, we may also write the dynamics directly for the SU(iV) matrix U and not in 
terms of the eigenvalues. After discretizing the Langevin time with stepsize e, this takes 
the form 

U(t + e) = exp [i\ a (eK a + y^fo)] U(t), (3.13) 

where A a (here a = 1, . . . ,N 2 — 1) are the traceless, hermitian Gell-Mann matrices, nor- 
malized as Tr A a A{, = 25 a \, and the drift terms are 

K a = -D a S v = i^TrXaU - if3 2 TvX a U~\ (3.14) 

Note that the complex process runs in the SL(iV, C) extension of SU(iV) and thus, e.g. in Eq. 
(3.13), U £ SL(iV, C) and K a G sl(iV, C). In this formulation there are iV 2 — 1 stochastic 
variables, instead of N — 1. In ref. [7] it was shown how the additional noise variables 
generate the correct drift for the remaining degrees of freedom when U is diagonalized 
after each update. This formulation was considered in refs. [6, 7, 17, 22]. 

3.3 Classical flow 

The role of the reduced Haar measure can be seen by studying the classical flow diagrams. 
Since the classical flow for the SU(3) spin model with angles cpi t 2 takes place in four real 
dimensions, it is hard to give a proper graphical representation. Instead we present as 
an illustration the classical flow for the SU(2) one-link model, in which the reduced Haar 
measure plays a similar stabilizing role. 

Explicitly, we compare the U(l) and SU(2) case, with 

Z m) = r dxe pcosx , Z sv(2) = [ W dx sin 2 xe pcosx , (3.15) 



U(l) " «« > ^SU(2) 

J —TV J —TV 

with complex (3 = /3r + if}\. The real and imaginary parts of the drift K, after complexifi- 
cation x — >■ x + iy, are 

2dsin(2x) 

K R = — Pr sm x cosh y + pi cos x smh y — — , (3.16) 

cos(2x) — cosn(2y) 

2dsinh(2y) 

K\ = — pr cos x sinh y + pi sin x cosh y H — — — -. (3-17) 

cos(2x) — cosh(2y) 

In the U(l) case, the contribution from the Haar measure is absent and d = 0; in SU(2), 
d= 1. 

Let us first consider the contribution from the Haar measure. It is easy to see that 
when \y\ is large, 

Kr — > 0, K\ = — 2sgn(y) (Haar measure only), (3.18) 

i.e. the contribution is purely restoring. We also note that at the origin and at (±7r, 0) 
the drift force from the Haar measure becomes singular; this invalidates potentially the 
formal argument for correctness of the complex Langevin approach method, because it 
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Figure 2. Classical flow diagrams: drift from the Haar measure only, in SU(2) (left); drift from 
the action only, in U(l) and SU(2), with j3 = 1 (right). 



might produce boundary terms at the singularity. Practically, however, this does not seem 
to cause problems, as the equilibrium measure appears to vanish sufficiently strongly at the 
origin (although it does require the use of an adaptive stepsize algorithm in simulations 
[23]). Finally, there are stable fixed points at (x,y) = (±7r/2,0). These findings are 
illustrated in Fig. 2 (left), where the drift from the Haar measure only is shown. 

Let us now consider the contribution from the action by itself, i.e. the drift in the 
U(l) model (d = 0), see Fig. 2 (right). This drift has an attractive fixed point at the 
origin and repulsive fixed points at (±7T, 0). Due to these repulsive fixed points, the real 
manifold (y = 0) is unstable against small fluctuations in the y direction. Therefore even 
with stochastic kicks in the x direction only (real noise), fluctuations in the imaginary 
direction, e.g. due to a small but nonzero f}\, will immediately lead to an exploration of 
the complexified space, resulting in trajectories that may take large excursions in the y 
direction, due to the character of the repulsive flow around x = ±7r. 

In the nonabelian case, one has to add the contribution of the Haar measure to that of 
the action. The resulting flow is shown in Fig. 3 (left) for real (3 = 1. The singular flow at 
the origin and at (±7r, 0) remains. We observe that the repulsive fixed points have moved 
away from the real manifold to (±7r, ±y*), while the attractive fixed points remain on the 
real manifold at (x* , 0) , where 




d? did 2 
cosx* = ~ + \jl + -j, coshy* = - + Wl + -j, (3.19) 



(for (3 = 1 this yields x* = ±1.14, = ±1.53). Note that the presence of the reduced 
Haar measure (d = 1) is essential. Since the repulsive fixed points have moved away, the 
real axis is stable in SU(2) for (3 values that are not too large. This can be seen by a linear 
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Figure 3. Classical flow diagrams for SU(2) for /5 = 1 (left) and /3 = 1 + 0.5i (right). 



stability analysis, which yields (for real f3) 

Ad 

y = -\ y a = /3cosx + - — — . (3.20) 

1 — cos(2x) 

This is stable (A > 0), as long as fi/d < 5.196. In fact, in the entire strip bounded by ±y*, 
the flow is directed towards the real manifold. 

Let us now see how a complex coupling changes the flow, see Fig. 3 (right) for /3 = 
1 + 0.5i. Even though the stable fixed points move away from y = [to ±(1.12, —0.158)] 
and the unstable ones from x = ±7r [to ±(2.79, —1.41)], we note that in the region around 
y ~ the flow is still directed to the real axis. This means that if initial conditions are 
chosen close to y = 0, the complex stochastic process will take place in a strip around 
y = 0. The probability distribution P(x, y) will be strictly zero outside this strip and the 
theoretical foundation of the complex Langevin method is justified [13, 14]. For larger 
(complex) f3 the stable fixed points will move out further into the complex plane and the 
dynamics will eventually no longer be confined to a strip. The stability analysis as carried 
out here is then no longer applicable. 

In the SU(3) case similar conclusions hold. The presence of the reduced Haar measure 
is essential in stabilizing the real manifold (for nearly real dynamics, i.e. when /3i,/x <C 1). 
Instead of showing four-dimensional flow diagrams, we look at the linear stability of the 
real manifold in the case that h = and (3 is real. We write (j) a — > (j^ + iy a and expand to 
linear order in y a (a = 1, 2). We find 



A = \ a a° • ( 3 - 21 ) 






The contributions from the action and the reduced Haar measure are denoted with U resp. 
H and we write A = A u ± A H . We find 

A u a = 2(3 [cos(^) + cos^f + $)] , A U Q = 2/3cos(0f + cf>f), (3.22) 
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Figure 4. SU(3) one-link model: (TrU k ) vs. 7 for k = 1,2,3, at fi eS = 2,/x = 1 and /i = 0.02 
(left) and 0.25 (right), using the two- and three-angle and the matrix formulations. At h = 0.02, 
agreement with the exact results is seen for all 7 in the two- and three-angle case; in the matrix 
case there is poor convergence at specific values of 7. At larger h = 0.25, there is poor or wrong 
convergence also in the two- and three-angle case. 



and 



A? = =■ CSC Z ( 2^—^ ) + 2csc z ( ZZ±_L±. ) + Z csc 



2 v 2 / V 2 / 2 V 2 

A? = _ 1 ^ (g-gj + csc2 ^2<+^ + csc2 ^< + 2^ (3 23) 

where b ^ a and csc z = 1/ sin 2. 

It is easy to see that the eigenvalues of by itself are positive for all values of 0^ 2 ' ^ ne 
real manifold is stable when only this contribution to the drift is considered. On the other 
hand, the sign of the eigenvalues of A u depends on 0^ 2 > leading to unstable regions. In 
the combined linearized drift, with A = A u + A H , inspection shows that the real manifold 
is stable for j3 < 1.75, and (linearly) unstable for larger f3. We note that this bound covers 
most of the parameter space indicated in Fig. 1, except deep in the ordered phase, where 
(3 can be larger (note that /3 here refers to (3 e g = 6/3\u\ in Fig. 1). Here we remark that the 
linear stability analysis is probably too restrictive, since no problems were encountered in 
this region. In fact, in the ordered phase complex Langevin dynamics was seen to perform 
very well [12]. 

3.4 Numerical results 

We have solved the complex Langevin equations numerically for a number of parameter 
values, using the two- and three-angle and the matrix formulations. Here we give a brief 
summary, illustrating the analytical findings discussed above. 

Firstly, we find that, although they imply different processes, the two- and three-angle 
formulations give identical results (within the statistical error) and agree with the exact 
results for small values of (3 and h, see e.g. Fig. 4 (left) for j3 = 2, h = 0.02, fi = 1. 
Since ju^O, the dynamics does not take place on the real manifold but in the complexified 
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Figure 5. Expectation values (Tr U k ) as a function of k for various values of /3 and \i in the SU(3) 
one-link model (left) and in the three-dimensional model on a 10 3 lattice. On the left, exact results 
are indicated with plus symbols; these can hardly be distinguished from the complex Langcvin 
results. 

configuration space. Nevertheless, the drift from the reduced Haar measure has a stabilizing 
effect, which constrains the dynamics. For larger values of the effective couplings /3i,2 (a*) 
(which can be achieved by increasing /3, h or fi), we find that the dynamics may break 
down for a range of 7 values. This is illustrated in Fig. 4 (right) for f3 = 2, h = 0.25, jjL = 1. 
A breakdown can be characterized by large fluctuations (e.g. around 7 ~ ±1.25) or by 
convergence to the wrong result (around 7 ~ ±7r). This happens at parameter values 
where the dynamics is not clearly dominated by a stable fixed point in the complexified 
field space. Interestingly, when the one-link model is viewed as an effective model for the 
three-dimensional SU(3) spin model (or full QCD), the combination of parameter values 
where the dynamics breaks down appears to be in a region of parameter space which is 
less relevant for those; recall the discussion around Fig. 1. Therefore this breakdown does 
not undermine the results obtained in the three-dimensional case. 

In the matrix formulation, we observe that the dynamics is more subtle and can break 
down earlier. This process is distinct from the previous ones in that a different complexified 
field space is explored, with apparently more possibilities for unstable trajectories. This is 
visible in Fig. 4 already at h = 0.02 for k = 3 at specific 7 values and more so at larger 
h = 0.25 (here the result at k = 3 is not shown). It is of course mandatory to control 
the dynamics in this formulation in the context of gauge theories and this can be partly 
achieved using gauge cooling [24]. 

The mathematical foundation for complex Langevin dynamics to converge to the cor- 
rect result is that the distribution in the complexified field space falls off rapidly. In ref. 
[12] it was found that in the three-dimensional SU(3) spin model the distribution in the 
imaginary direction c/) 1 drops exponentially, P{$) ~ e - ^ 1 ', with b ~ 35 — 45. This left 
open the question what happens to expectation values of the form (TiU k ) with k large. 
These observables contain terms of the form e^ k< ^ cos(/c<^ R ) and the presence of the rapidly 
oscillating cosines should be taken into consideration. This is demonstrated in Fig. 5 where 
(Tr U k ) is shown as a function of k in the one-link model (left) and in the three-dimensional 



- 11 - 



model on a 10 3 lattice (right). In the one-link model, exact results are shown as well; agree- 
ment between the Langevin and exact results is observed. We conclude therefore that the 
higher moments decrease rapidly and that the falloff of the distribution in the imaginary 
direction is sufficient to achieve this. 

Our conclusion from the numerical experiments with two and three angles is that the 
Langevin dynamics is controlled and yields the correct results when the contribution from 
the reduced Haar measure is not overpowered by the contribution from the action itself. 
This is achieved by taking the effective couplings not too large. For larger (3 or h values 
however, the real manifold is no longer stable and this implies problems for the complex 
Langevin process. Fortunately, this happens for parameter values in the effective one-link 
model, which do not seem to be relevant for the three-dimensional (or four-dimensional) 
original model. 

Above the one-link model was used as an effective model for the SU(3) spin model. 
However, we can also see it as an effective model for lattice QCD in the heavy dense 
approximation, cf. refs. [20-22, 24]. In this case the effective coupling f3 represents the 
contributions from the staples attached to a link, and instead of Eq. (2.3) the heavy dense 
determinant itself appears in the three-dimensional effective action: 



with h = (2k) * and V* the (conjugate) Polyakov loops. Using the determinant in 
the corresponding one-link model appears to have a stabilizing effect and the two- and 
three-angles simulations are even better behaved. 

4 Stabilization via generalizations of the CLE 

A general possible route for stabilization uses the fact that there is a wide range of pos- 
sibilities to modify the complex Langevin equation (CLE). Such a modification, if chosen 
appropriately, can improve the falloff of the equilibrium distribution, which was identified 
in refs. [13, 14] as the essential prerequisite for a correct CLE process. Some of those 
modifications, like kernels, were known already in the 1980s [25-27], others such as vari- 
able transformations, were encountered just recently. In fact the reduction to the Cartan 
subgroup discussed in the previous section falls into this general category. In the present 
section we analyze the possible modifications more generally and find an equivalence be- 
tween the approaches using a kernel and coordinate transformations. 

The crucial point is that we apply the CLE only for averaging holomorphic functions. 
Different distributions in the complexified configuration space can give the same expecta- 
tion values for those holomorphic functions; this is true for the probability density evolving 
under the Fokker-Planck equation as well as the stationary equilibrum distribution. In Ap- 
pendix A we give general classes of processes and equilibrium measures which are equivalent 
in this sense. In addition there is the well-known freedom of using different processes having 
the same equilibrium measure. 




X 



(3.24) 
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Let us first describe the use of coordinate transformations. We have to distinguish two 
kinds of transformations: 



(a) transforming the CLE process, i.e. the description of the trajectories; 

(b) transforming the variables in the functional integral. 

The first possibility does not change anything essential, since it is only a transformation 
of the description; the second choice, however, leads to a modified process due to the 
appearance of a Jacobian. It turns out that after changing the description of this modified 
process back to the original coordinates one has in fact introduced a kernel. 

In one variable this works as follows: let 

x = x(u) (4.1) 

be an invertible smooth variable transformation. To avoid confusion we rename the action 
S(x) when considered as function of u by 

S(u) = S(x(u)) . (4.2) 

We have 

Z= dxe~ s ^= / due~^ u \ (4.3) 

with 

S cS (u) = S(u) - In J(u), J(u) = (4.4) 

The new drift becomes 

K{u) = -S' cS (u) = -S'(u) + J'(u)/J(u), 

and the new CLE 

u = -S'{u) + j'(u)/J(u)+7](t). (4.5) 

Let us first study the stability of the real manifold under small imaginary fluctuations in 
the original formulation. We write x — >■ z = x + iy and linearize in y to find 

y = -X(x)y, A(z) = S"(x). (4.6) 

Only if X(x) > for all x is the real axis stable against complex perturbations. 

By choosing a clever transformation one can stabilize a situation which is unstable due 
to a nonpositive X(x) or a repulsive fixed point. After the change of variables to u and 
writing u — > u + iv , one finds, to first order in v , 

v = -\(u)v, (4.7) 

where now 

A(n) = SM = S"(u) - J"(u)/J(u) + [J'(u)/J(u)] 2 . (4.8) 
The contribution from the Jacobian may now stabilize the real manifold. 
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Figure 6. Classical flow diagrams in the Gaussian model with a = — 1+i in the original formulation 
(left) and after the coordinate transformation x — it 3 (right). 



A striking example is given by the Gaussian model, with the action 



S 



--ax 



a = a + ib, 



(4.9) 



which is well defined as long as Recr = a > 0. The origin is an attractive fixed point in this 
case. But for a < not only does the integral become ill-defined, after complexification 
also the fixed point at the origin becomes repulsive, see Fig. 6 (left). 

A simple transformation x(u) = u n (n odd) may remedy this situation, as shown 
in Fig. 6 (right) for a = + i and n = 3: while the additional term in the drift, 
J'(u)/J(u) = (n — yields a singularity at the origin, it also partly stabilizes the real 

axis, since 



A(u) = an(2n - l)u 2 ^ n - l) 



+ 



n — 1 



(4.10) 



and the second term is always positive. Most importantly, however, the change of vari- 
ables leads to the appearance of stable fixed points in the complex plane. We find that a 
numerical solution of the CLE in u converges and yields for (x 2 ) = (u 6 ) 



1 a + ib 



(4.11) 



a a 2 + b 2 ' 

even for a < 0, see Fig. 7 (left). This answer is correct in the sense that it is the analytic 
continuation of the result valid for a > 0. 

The same transformation also allows to evaluate the quartic 'Minkowski' integral [28, 

29] 

Z = r dx exp (-^x 4 ) , (4.12) 



defined by analytic continuation, giving 



^ ' Vaf(|) 1 h 



(4.13) 
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Figure 7. Complex Langevin results for the correlator (x 2 ) in the Gaussian model with a = a + i as 
a function of a (left) and in the quartic model as a function of A (right), both after the transformation 
x = u 3 . The lines indicate the expected, analytically continued results. 



see Fig. 7 (right). 

As the third example, we take a more realistic case, encountered in the previous section, 
namely 

Z= r dxe^ cosx 03 GC). (4.14) 



We would like to stabilize this by a variable transformation that casts the integral in the 
form of an SU(2)-like integral, see Eq. (3.15). This can be achieved by the change of 
variables, 

x(u) = u — — sin(2n), — ir < u < tt, (4-15) 

since the Jacobian is 

j( u ) = ^M = 2 sin 2 u. (4.16) 



du 

Simulations in the new formulation yield results shown in Fig. 8. The results are in line with 
the findings from above. For small j3 = 0.3 (left), the Jacobian indeed has a stabilizing effect 
and the exact results are reproduced. For larger /3 = 1 (right), the unstable contribution 
from the action takes over and exact results are only reproduced for 7 ~ 0, it. In the 
three-dimensional XY model, a variable change along these lines gave only partial success, 
unfortunately. 

We now want to show that these transformations are equivalent to the introduction of 
a kernel. The CLE with a kernel can be written as 



z = H 2 (z)K(z) + 2H(z)H'(z) + H(z)r](t), 



(4.17) 



with z = x + iy and a holomorphic function H (see refs. [25-27], where H 2 is called the 
kernel). Separating real and imaginary parts, Eq. (4.17) becomes 



x = K x + ReHr](t), 



i, = K y + Im Hrj(t), 



with 



K 



-H 2 X7 Z S + V Z H 



2 



K x = ReK 



K y = ImK. 



(4.18) 
(4.19) 
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Figure 8. Complex Langevin results for the correlator (cosx) in the U(l) model, with complex 
coupling /3e 47 as a function of 7 at f3 = 0.3 (left) and f3 = 1.0 (right), after the transformation 
x = u — i sin(2w). 



Let us now change the description of the transformed process back to the original variables. 
This requires the use of Ito calculus and the chain rule; one obtains the process 

+ (uY 1 V (t), (4.20) 

This can be recognized as the original process modified by a kernel using H = 1/u' = J. 
The fact that H strictly speaking is not a holomorphic kernel because of its singularity 
at the origin is of no consequence in the cases discussed before, because the equilibrium 
distribution vanishes at this point. 

The reduction of the SU(2) and SU(3) effective models to the Cartan subgroup dis- 
cussed in Section 3 is an example of singular variable transformation; the same is true 
about the relation between the iV angle and the (N — 1) angle formalisms discussed there. 
In Appendix A these issues are discussed in a more systematic way. 

5 Conclusion 

Complex Langevin dynamics can handle theories with a complex action, even when the 
sign problem is severe, but success is not guaranteed. In this paper we contrasted the 
apparent success in the three-dimensional SU(3) spin model, an effective model for QCD 
at nonzero chemical potential, with the observed failure in the three-dimensional XY model 
in the disordered phase, by constructing and analyzing effective one-link models. 

We found that a crucial role is played by the presence of the nontrivial Haar measure 
in the nonabelian case. The contribution from the reduced Haar measure to the drift 
constrains the dynamics in the complexified field space: it is purely restoring and directed 
towards attractive fixed points on the real manifold. The presence of a singularity in the 
drift at the origin does not cause problems. On the other hand, the drift from the action 
yields repulsive fixed points, which can be located on or close to the real manifold and are 
responsible for large excursions in the complexified field space. As we have shown earlier, 



(u') 2 



S'(z) 



2u" 
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an uncontrolled exploration of the enlarged field space results in insufficient falloff of the 
equilibrium distribution, which causes the formal justification of the complex Langevin 
equation to break down and also leads to wrong results in practice. 

In the nonabelian case, the contributions from the reduced Haar measure and action 
balance each other, resulting in a controlled Langevin process. We showed that the real 
manifold is stable against small complex fluctuations. This ensures analyticity of the 
observables in e.g. /i 2 . By increasing the importance of the contribution from the action, 
the Langevin process can be made to break down. However, the parameter values in the 
effective one-link models for which this happens do not seem to be relevant for the three- 
dimensional SU(3) spin model. The results of this paper provide further justification for 
the conclusion drawn in ref. [12], in which complex Langevin was argued to be applicable 
by carrying out a series of assessments, without relying on knowledge of the 'exact' result. 

In the (abelian) XY model the Haar measure is trivial: hence there is no restoring 
component and the real manifold is unstable against small complex fluctuations. Since 
these fluctuations are introduced as soon as fj, 7^ 0, simulations when fj, — >■ and /i = 
do not agree. This is indeed what has been observed in the disordered phase of three- 
dimensional XY model and is now fully explained. Deep in the ordered phase, the theory 
becomes effectively real and there is no problem for complex Langevin dynamics. 

The important lesson is that it is desirable for the real manifold to be stable against 
complex fluctuations in order for complex Langevin dynamics to work (even though strictly 
speaking it is neither sufficient nor necessary). This provides a useful benchmark to address, 
given that there is a wide range of options to modify the complex Langevin approach 
(some of them were discussed here). This is especially relevant for gauge theories, where 
the complexification leads to an extension from SU(iV) to SL(iV, C). Here it is possible to 
stabilize the real group manifold via gauge cooling, at least partly [24]. 
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A Generalizations of Langevin dynamics 

In this appendix we investigate systematically possible generalizations of complex Langevin 
(CL) dynamics that are at least formally equivalent to the original one. For the formal 
equivalence (i.e. ignoring problems of slow falloff and possible boundary terms, see refs. 
[13, 14]) it is sufficient to verify that the (complex) Fokker-Planck (FP) operator has 
exp(— S) as a zero mode. In refs. [13, 14] it is also explained how the equilibrium distribution 
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of the real FP operator acting on functions of the complexified configuration space becomes, 
when evaluated on holomorphic observables, equivalent to that zero mode. 



A.l Nonuniqueness of measures for holomorphic observables 

For simplicity we consider the 'flat' case of R . The goal of the complex Langevin equation 
(CLE) is to find a probability measure on C , 

P(x,y)dxdy x,y£R N , (A.l) 

which for holomorphic observables O is equivalent to a complex measure on M. N , 

p(x)dx, (A. 2) 

in the sense that for all O, 

dxdyP(x,y)0(x + iy) = dxp(x)0{x). (A.3) 



The question of the existence of such a measure P was already discussed by Weingarten 
[30], but here we are interested in the nonuniqueness. 

Clearly P is not uniquely determined by the requirement (A.3). First one has to 
recognize that a precise answer to this nonuniqueness depends on specifying the space of 
holomorphic observables for which Eq. (A.3) is supposed to hold. We consider spaces T~Lf 
defined in terms of some growth condition 

sup a{x, y)\0(x, y)\ < C < oo, (A.4) 

for some positive weight function a. Examples for / are / = exp((rr 2 + y 2 ) a ), with x 2 = 

X~^ ~\~ • • • X . CtiC. 

Obviously Pi and P2 are equivalent in the sense above if and only if for all O E Hf 
Q = P\ — P2 satifies 

/ dxdyQ{x,y)0(x + iy) = 0. (A. 5) 

Jc N 

In other words, Q has to be orthogonal to the intersections of the kernels of the Cauchy- 
Riemann operators 

(';, <),-, ■ i<) !h . (A.6) 

Ignoring possible convergence problems or boundary terms, Eq. (A. 5) seems to say only 
that Qdxdy is a signed measure satisfying 

N 

Q{x,y) = Y J C 3 H 3 {x,y), (A.7) 
i=i 

where the Hj are some complex- valued functions (of bounded variation), because rolling 
the CR operators Cj over to act on the observables would yield zero. Those conditions are 
satined if there are real functions Gj(x,y) satisfying 

N 

Q(x,y)=^2A j G j (x,y) (A.8) 
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with Aj = d 2 x . + d 2 ... 

The convergence problems have to be taken seriously; we are only sure that they are 
absent if QO is absolutely integrable. Depending on the class of holomorphic observables 
chosen, one has to demand corresponding decay properties of the Gj. 

To make use of Eq. (A. 7) for a cure of convergence to the wrong limit, one would need 
to do the following: first find a Q that is such that it improves the decay of P + Q at 
infinity (this is of course only possible if Q shows that same rate of decay as P), then find 
a modified CL process having P + Q as its equilibrium measure. 

A. 2 Neutral modification of the CL process 

Next we ask how one can modify the CL process without changing the evolution of holo- 
morphic observables under the Langevin evolution; we call this a neutral modification. If 
the Langevin operator is given by 

L=[V x + K x ]V x + K y Vy, (A.9) 

we want to see if we can replace it by an operator of the general form 

3 3 

with coefficient functions Fj,Gj,R x ,R y (F,G > 0) in such a way that when applied to 
holomorphic observables O, L' becomes equal to L, because d Vj O = id Xj O (j = 1, . . . N). 

It is not hard to see that this requires (1) Gj = Ff (j = 1, . . . A); (2) R x = R y = 0. 
So the modifying operator is of the form 

I.,„ »/(^. • (A.ll) 

j 

corresponding to the modified CL process 

±j = K Xj + v 1 • F'f Vj (t) , /), = K Vj + Fj m (t) . (A. 12) 

A well-known special case arises by choosing 

F J (x,y) = N I , j = l,...N, (A.13) 

which we refer to as 'complex diffusion'. Here we have discovered a generalization: 'configuration- 
dependent complex diffusion'. But it is not clear how much this can help in solving the 
problems of the CLE, in view of the fact that it was found that introducing 'imaginary 
diffusion' (or complex noise) actually tends to create problems by leading to slow decay in 
the imaginary directions and convergence to the wrong limit [13, 14]. So Fj should at least 
be chosen in such a way that it vanishes sufficiently fast for \yj\ — > oo, y = 1, . . . , N. 
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A. 3 Nonneutral modifications of the CLE: holomorphic (matrix) kernels 

A way of modifying the CL process that was discussed already in the 1980s is the introduc- 
tion of a so-called kernel [25-27] . It is nonneutral in the sense that it does not preserve the 
Langevin evolution of holomorphic observables (but is supposed to preserve the long-time 
averages of holomorphic observables). For one dimension it amounts to the following: 

Let H be holomorphic on C (H corresponds to what is called \f~K in the older litera- 
ture). We denote 

ReH = R, lmH = I. (A.14) 

The corresponding CLE becomes 

x = K x + Rri(t), y = K y + Ir}{t), (A.15) 

where 

K = -H 2 S' + 2HH', K x = ReK, K y = lmk. (A.16) 

It is easy to find choices for the kernel that eliminate the drift force altogether or change its 
sign, but the resulting process suffer from steeply rising diffusion coefficients and dubious 
stability. 

A special case that does not suffer from these problems is a constant (complex) kernel 
H, with ReH 2 > 0. In this case the kernel just amounts to multiplying the time coordinate 
by a factor H 2 and there are no consistency problems. This simple device has been shown to 
improve the situation in certain cases with quadratic action in refs.[26, 27]; the stabilization 
occurring there is actually equivalent to the one discussed in Section 4. 

For more than one variable there is a more general way of introducing kernels which 
we call matrix kernels. In short, take a holomorphic Hj k (x\ + iy±, ■ ■ ■ ,x n + iy n ), then 

®j = K x ,j + Re ^ H jk rj k (t), yj = k yJ + Im ^ H jk t] k (t), 

k k 

where 

&j = E {-{H T H) kj {V Zk S) + V Zk {H T H) jk ) , 

k 

k x j = Re Kj , k yJ = Im Kj . (A. 17) 

To see that this is (formally) correct, note that the complex Fokker-Planck operator is 
l h = E {Vj(H T H) kj [V fc + (V k S)}} . (A.18) 

k,j 

It is manifest that exp(— S) is a (hopefully unique) zero mode of L T . 

We remark that we can actually combine the two previous modifications (kernel and 
imaginary diffusion). For simplicity we only describe the case of one dimension: let again 
H be a holomorphic function and F a real valued function on A4 C . Then the modified 
Langevin operator 

L H = (OR 2 + F 2 )V X + it) V, + ((/ 2 + F 2 )V y + K y ^ V y + 2RIV x V y (A.19) 
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reduces on holomorphic observables to 



L H , F = H 2 A Z + KV Z , (A.20) 

as before. 

We have to find a process corresponding to these operators. This problem does not 
have a unique solution, but the following stochastic differential equation turns out to give 
the right Lh,f- 

. ~ Ry/ff + P + F* FI 

X = K x H ; Vl(t) 1?2 t , 

IV R 2 + P + F 2 FR . A , 

y = K v + ; Th(t) + —= m(t), A.21) 

y VR 2 + 1 2 ' VR 2 + I 2 

where 771, 772 are two independent white noises. Note that these equations reduce to the 
previous ones for the special case F = 0; for F = Nj and R = 1, I = one recovers the 
standard complex noise. 

A. 4 Nonneutral modifications of the CLE: Coordinate transformations 

As remarked in Section 4, coordinate transformations are closely related to kernels. In 
fact they lead to a subset of the general matrix kernels discussed above, but since they are 
sometimes easier to handle we discuss them in some detail. To avoid heavy notation, we 
treat real actions and comment about the complex case in Subsection A. 4. 6. 

A. 4.1 General remarks on coordinate transformations 

As stated in Section 4, we have to distinguish two types of coordinate transformations: 

(I) Transforming the descriptions of the trajectories; 

(II) Transforming the functional integral. 

Choice (I) does not change the process, only its description, but choice (II) is nontrivial. 
The difference can be seen clearly in a simple linear one-dimensional example: let 

u = ax {a / 0) (A. 22) 

and consider the real Langevin process corresponding to the action S{x) = S(u/a) = T(u): 

x = -S'(x) +v(t). (A.23) 

Then (I) means 

(I) : ii= -a 2 T'{u) + arj(t), (A.24) 

whereas (II) means replacing the measure proportional to exp(— S(x))dx by the trans- 
formed one proportional to exp(—T(u))du and setting up a Langevin process that has the 
transformed measure as its equilibrium measure. We find 

(II): u = -T'(u) + r](t). (A.25) 

So in this simple example the difference between (I) and (II) is just a different time scale. 
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A. 4. 2 Nonsingular (invertible) linear transformations 

We now consider a less trivial case of a Langevin process in R , using vector notation: 

x = -VS(x) + ?](t), (A.26) 
and the transformation by a nonsingular matrix A 

u = Ax, r(«) = 5(A _1 u). (A.27) 

Scheme (I) yields, simply by acting with A on both sides of the LE, 

u = -AV x S(x) + Adr] for it = Aa;. (A.28) 
To re-express the drift in terms of T(u), we use the chain rule 

V x T(Ax) = A T V u T(u) for u = Ax, (A.29) 

which leads to 

(I) : it = -AA T VT(u) + Arj(t), (A.30) 

whereas (II) gives 

(II) : u = -VT(u) + r](t). (A.31) 

The two processes now no longer differ just by different times scales. It is also instructive 
to compare the Fokker-Planck operators for the two processes. In matrix notation 

(I) : L T = VAA T (V + (VT)) (A.32) 

and 

(II) : L T = V(V + (VT)). (A.33) 

It obvious that both operators have the same equilibrium solution. The transition from 
scheme (II) to scheme (I) can be viewed as the introduction of a constant 'matrix kernel' 
A; in this sense it is more general than the standard introduction of scalar kernels. 

A. 4. 3 Singular (noninvertible) linear transformations 

If the action has a certain symmetry, noninvertible linear transformations can be used to 
reduce the number of variables to a subset; conversely one can use them to introduce a 
redundant set of variables, which is sometimes useful. 

Here we show explicitly how this can be used to justify the process on the Cartan 
subgroup of SU(N), using N instead of N — 1 angles, as discussed in sec. 3.2. We have 
to consider singular linear transformations A that map M. N into a lower dimensional space 
R M (with M <N): 

u = Ax. (A.34) 

A has a null space V Q C M> N and an orthogonal complement V± . We assume that the action 
is invariant under shifts in the V$ directions, i.e. 

S(x + u) = S(x) whenever u G Vq. (A. 35) 
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Then the drift VS will leave V± invariant. 

Specializing now to the case of SU(3) and going from the 3 angle to the 2 angle 
formulation, we may think of factorizing A as 

A = QB, (A.36) 

where B is a singular symmetric N x N matrix annihilating Vq and mapping V± onto itself. 
Q is then a one to one mapping of Vq onto M M . 

In the case of Eq. (3.10), with x resp. u in Eq. (A. 34) denoted with z resp. cj) in sec. 
3.2, B maps the z's onto the subspace with <j>\ + fa + 4>3 = in terms of the 3 angles 
01 ) 4>2 , 03 j so we have 

/ 2 -1 -l\ 
-1 2-1 , (A.37) 

V-i-i 2 ; 

and Q maps onto the two independent angles <f>\,4>2i he. 



B = i 
3 



so finally 

A = OR = 

3 V — 1 2-1 



,4 = QS = - ( 2 1 1 ) . (A.39) 



The LE for the variables x £ is 

i = -VS(i) + j/(t), (A.40) 

where r\ consists of independent white noises, hence does not leave the subspace Vq 
invariant. The remainder of the analysis is almost identical to the one for nonsingular A: 
the LE for u = Ax (scheme (I)) can at first be written as 

u = -AV x S(x) + Arj(t), (A.41) 

but we have to re-express the drift —'VS(x) in terms of the transformed variables u. First 
we define 

T(u) = S(x) whenever u = Ax, (A. 42) 

which is well defined because of the invariance of S(x). By the chain rule again 

V x T(Ax) = A T V u T(u), (A.43) 

so the LE for u becomes 

(I) : u = -AA T V u T(u) + Arj{t) . (A.44) 
In scheme (II), on the other hand, the LE for u is simply 

(II) : u = V u T{u) + V (t) , (A.45) 
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where rj(t) consists now of M independent white noises. To see to which extent the two 
schemes are equivalent, let us look at the corresponding FP operators: 



(I) : Lf = V T U AA T {V U + (V U T)) (A.46) 

and 

(II) : Ll = V U (V W + (V„r)) . (A.47) 

Again it is obvious that both schemes lead again to the same equilibrium measures, but 
clearly the processes are quite different. 

A. 4. 4 Invertible nonlinear coordinate transformations 

The treatment of nonlinear coordinate transformations is quite similar to the linear case, 
except that one has to use Ito's calculus (see for instance ref. [31]) in some places. Again 
we have to consider both schemes (I) and (II). We first consider invertible transformations 
of ~S, N and comment in the next subsection about singular (noninvertible) ones. 

Let / : W N —7- M. N be an invertible and continuously differentiable coordinate transfor- 
mation. If (in vector notation) 

u = f{x) (A.48) 
and x(t) evolves according to the Langevin equation 

x = -VS(x) + ?](t), (A.49) 

we define again 

T(u) = S(f-\u)). (A.50) 
Then by the Ito calculus, u evolves according to the LE 

. {k^ + H + ^ (t) /- (A - 51) 

To express this LE in terms of u and T(u) we again use the chain rule 

dSju) _ydu r dT(u( X )) ^ 
dXj Q x . Q Ur 

If we define the matrix A by 



A i3 - (A.53) 



the LE for u becomes in matrix notation 

(I) : u = -AA T V u T{u) + A x u + Ar](t), (A.54) 

which clearly reduces to Eq. (A. 44) in the linear case. 

In scheme (II) the transformation of the functional integral involves the Jacobian 

J=^4 = detA" 1 , (A.55) 
o(u) 
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so the action T(v) receives an extra term 

Tj = - In J = bidet A (A. 56) 

With 

T cS = T + Tj = T + ln det A, (A.57) 
we get the Langevin equation for scheme (II) 

(II): u = -V u T cS (u) + rj(t), (A.58) 

with the rj(t) denoting again N independent white noises. 

To see that the two schemes are equivalent as far as their equilibrium measures are 
concerned, we again compare the Langevin and FP operators of the two schemes. We first 
note a general fact (see for instance ref. [31]): if 

x = K(x) + A(x)rj(t), (A.59) 

then by the Ito calculus 

L = Y / (AA T ) lk (x)d i d k + K(x)-V (A.60) 

Lk 

and 

L T = Y,didk{AA T ) lk (x) - V • K(x). (AM) 

i.k 

Applying this to our two Langevin equations above we find for scheme (I) 
and for scheme (II) 

(II) : Ll = Vl [V u + (V u T cS (u))} . (A.63) 
We claim that the FP operator in scheme (I) can be rewritten as 

(I) : Lf = V T U AA T [V u + (lV u T cS (u))} , (A.64) 

again making the equivalence of the two schemes obvious. To prove our claim we only have 
to show that 

^(AA^^- In det A = -A xUi + ^ JL(AA T ) ih ; (A.65) 

k ^ k ^ 

this follows again using the chain rule: 

T ^-(AA T ) ik = y p- A (^P). (A.66) 
z —' ouk fr 1 out ox\ \ dx r ox r ) 

The last expression is equal to 

c 

duk dx\dx T dx r ' duk dx r dxidx 



y J dxi d 2 U t du k | dXj dUj d 2 U k ■ 



k.l.r 
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here the first term simplifies to A x Ui, whereas the second term equals 

£(X4 r ) ifc ^-bidetA, (A.68) 

k % 

so Eq. (A. 65) is correct. 

As in the linear case the transition from scheme (II) to scheme (I) can be viewed as 
the introduction of a 'matrix kernel' A (see Eq. A. 53), which is a generalization of the 
standard introduction of scalar kernels. Again we can transform scheme (II) back to the 
original coordinates using scheme (I) and obtain a process with a matrix kernel A~ l . 

But note that the matrix kernels obtained here are not as general as before, since the 
matrix = dxi/duj has to satisfy an integrability condition, which was not required of 
the matrix H in sec. A. 3. 

A. 4. 5 Singular (noninvertible) nonlinear transformations 

As in the case of singular linear transformations, in this case we have to require some 
additional symmetries of the action. A simple example is the action 

S = -p\xU, (A.69) 

with U an element of SU(iV). The symmetry in question is the invariance of S under 
similarity transformations 

U ' — y V~ X UV, (A.70) 

which allows the transition to the Cartan subgroup of SU(iV). The map from the N 2 — 
1 parameters of SU(iV) to the N — I Cartan angles is clearly a singular noninvertible 
transformation. 

More generally, in a lattice gauge model gauge invariance allows to reduce the process 
to one remaining in a lower-dimensional gauge-fixed configuration space. 

A. 4. 6 Remark on complex actions 

If S(x) is complex, at least on the formal level (ignoring possible problems with boundary 
terms under integration by parts) everything goes through as before. 

The formal proof of correctness given in refs. [13, 14] only uses the fact that the complex 
FP operator annihilates the complex measure exp(— S)dx. This is manifestly the case for 
the most general case of a holomorphic matrix kernel, see Eq. (A. 18). 

Here in principle we have to limit ourselves to transformations that are bi-holomorphic 
maps from z = x+iy tow = u+iv. In practice it seems that this restriction is too severe and 
certain singularities can be tolerated - provided we are lucky and the equilibrium density 
vanishes at those points. The formal proof of equivalence again only uses the equivalence 
of the complex Fokker-Planck operators - the Langevin evolutions in the space are again 
different in the two schemes. 
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